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Abstract 

m ■ 

^ ' We calculate the low-lying spectrum of new particles expected in a minimal 

op ■ model wherein supersymmetry breaking at ^ 100 TeV is indirectly transmitted 

to the Standard Model. We calculate the couplings of these particles relevant 
l/^ ' to their most nearly accessible experimental signatures, and estimate those 

• signatures quantitatively. Running of various couplings plays a crucial role in 

the phenomenology, specifically in generating an adequate Ni — Ir splitting. 
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I . 1 Introduction 

D , 

^ ■ Since the widespread realization that successful unification of gauge couplings can be 

achieved in models including the minimal supersymmetric extension of the standard 
^ \ model (perhaps augmented with complete SU{b) multiplets) anticipation of the 

5^ \ imminent discovery of direct evidence for real - as opposed to virtual - supersym- 

metry has become almost palpable. On the other hand, there has been no consensus 
regarding the form this discovery might take. Theoretical control over the many new 
parameters that arise even in minimal supersymmetric extensions of the standard 
model is not adequate to guide one through a bewildering number of choices that 
substantially affect the predictions for observable phenomena. 

Recently a class of models that potentially enjoy great predictive power has been 
identified . The leading idea of these models is that the fundamental breaking of su- 
persymmetry occurs in a "hidden sector" and is characterized by a scale A < 100 TeV. 
This breaking is supposed to be conveyed to our observable sector, including the mini- 
mal supersymmetric standard model, by a specific, highly symmetric interaction term 
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(to be discussed more precisely below). This is to be contrasted with the situation in 
a class of models, commonly referred to as supergravity models, that has been much 
more popular. In supergravity models the fundamental symmetry breaking scale is 
A ~ 10^ TeV, and the hidden sector communicates with the observable sector by 
gravitational-strength interactions. A decisive phenomenological difference between 
these classes of models arises because of the different role played by the gravitino 
G, which is essentially the Nambu-Goldstone fermion of spontaneous supersymmetry 

A 2 

breaking. Its mass is generically given as Mg = -p^^ — -, whereas its couplings to ordi- 
nary matter scale as ^ 0. In most supergravity models the gravitino is massive and 
very weakly coupled, and although it could be significant cosmologically, it plays no 
very direct role in ordinary phenomenology. In the (relatively) low energy messenger 
models the gravitino is exceedingly light, and although its couplings are very feeble 
they can be such as to cause the next lightest, photino-like particle A^^i to decay via 
Ni ^ 'J + G within a distance that might be difficult to resolve experimentally [Q. 
Generically, then, in these models one will find that pair-production of supersymmet- 
ric (R-odd) particles will lead eventually, perhaps after a cascade of decays, to final 
states containing 77 and missing energy. The recent observation of a dramatic 
event with a final state containing hard e"'"e~77 and missing transverse energy by the 
CDF collaboration therefore lends special interest to these models. 

In this paper we shall calculate the low-lying spectrum and interactions in a min- 
imal messenger model of the kind mentioned, in sufficient detail to estimate the most 
nearly accessible processes. We find a remarkably tight and specific pattern of con- 
sequences. In the course of the analysis we shall find that some apparent difficulties 
of the model are resolved by careful attention to the running of various couplings. 
Specifically, the tiny Ni — cr mass splitting one finds using bare couphngs, and which 
if correct would render the interpretation of the observed e"'"e~77 event problem- 
atic 0, is found to be greatly enhanced by their running (and by the inclusion of 
D-term effects). Similarly a possible difficulty with the apparent smallness of the 
Bii parameter, which if valid would imply the existence of a weak scale axion, is 
automatically repaired. 



2 Spectrum 

In the minimal messenger model (MMM) supersymmetry breaking is conveyed to the 
standard model particles by the matter fields (g', £') -|- (g', J ) having the quantum 
numbers 5-1-5 under SU{b). These fields are assumed to acquire supersymmetry 
breaking masses through their couplings to a standard model singlet field X which 
develops vacuum expectation values along its scalar as well as its F-components. 
Explicit models where this occurs dynamically have been constructed [Q. Typically 
these models contain the superpotential couplings 

W = \qqX + Xil'l'X (1) 
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which renders a non-supersymmetric spectrum for the (g', £') fields through Fx 7^ 0. 
The (g', £') fields, through their standard model gauge interactions, induce SUSY 
breaking masses for the R-odd scalars {q,t) and gauginos {B,W,g). The gaugino 
mases arise through one-loop diagrams and are given (in the approximation Fx ^ 

\i {Xf) as 



M,(A) 



a^(A) 
47r 



A 



Here A = Fx/ {X) is the effective SUSY breaking scale, which is close to the mass 
scale of the 5 + 5 fields. The scalar masses arise through two-loop diagrams and are 
given as 



m2(A) = 2A2 
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We have adopted the usual SU{^) normalization of the hypercharge coupling, ai = 
I ay and defined F as Q = T3 + In Eq. (H), C3 = 4/3 for colored scalars and zero 
for sleptons, C2 = 3/4 for SU{2) doublets and zero for singlets. 

The mass relations of Eqs. (@)-® will arise in a large class of models which 
might differ in the details of supersymmetry breaking; it is sufficient that SUSY 
breaking is conveyed to the standard model particles through the 5 + 5 fields as in 
Eq. (|I|). Thus the mass relations would hold even if SUSY is broken spontaneously 
by the O'Raifearteigh mechanism (F-type breaking). The 5 + 5 fields is of course 
the simplest extension of the standard model sector that preserves the successful 
unification of the gauge couplings. 

The relations in Eqs. (|D-(@) imply that there is a clean separation of scales among 
the SUSY particles, as their masses scale with the interaction strength. Thus the 
squarks are the heaviest, followed by the gluino, left-handed sleptons W, right- 
handed sleptons (^ij), and the 5-ino, which is the next-to-lightest R-odd particle 
(the lightest is the gravitino G). Their masses (neglecting RGE-running) are in the 
ratios 11.6 : 7.0 : 2.5 : 2 : 1.1 : 1.0. From an immediate experimental perspective 
the most interesting particles are the lightest ones, in, the chargino (which we shall 
denote by W) and the two neutralinos (iVi and A^2)- The latter are essentially pure 
gaugino eigenstates since their mixing with the Higgsinos are suppressed by the large 
fi parameter required in the model. 

The relations in Eqs. (0)-(^ receive significant corrections from the renormaliza- 
tion group evolution. This evolution has several interesting consequences: 

• It drives the mass-squared of the MSSM Higgs doublet that couples to the 
top quark, which is positive at the scale A, to a negative value near Mz and thus 
facilitates electroweak symmetry breaking. The mass-splitting between Hu and H^, is 
given by 

(4) 



2, 
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Figure 1: (a) Boundaries of interesting regions in the (M2,tan/5) plane: the region 
above the dark, sohd hne does not have a QED-preserving potential minimum; in the 
region above the light solid line, the tr is the NLSP; in the region to the right of the 
dashed line, the^ecay^ — > WNi is kinematically allowed; to the right of the dotted 
line the decay N2 — > Nih^ is kinematically allowed, (b) Contours of mj — m^^ in 
GeV, for /x < 0. 

Observe that Hu receives a mass correction proportional to m~, which is much larger 
than m?r . 

J^u 

• It lowers the masses of the {B, W), and raises the masses of (as well as £l 
and q), from their bare values at A. Consequently, the mass splitting between B 
and in is enhanced at the weak scale. We find that a mass splitting of 10 GeV at 
A typically becomes (20 — 25) GeV at Mz- Upon including, as well, the positive 
D-term correction to the mass of i^, the mass splitting becomes about 35 GeV for 
M2 — 200 GeV. (For definitiveness we have taken A = 100 GeV for the numbers given 
here.) Contours of this splitting are shown in Figure This difference has significant 
impact on the interpretation of the observed 6677 event at CDF. 

• It supplies a sufficiently large mass to the would-be axion of the model. The soft 
bilinear mass term Bfi for the Higgs field, {BuHuHd + H.c.) in the Higgs potential, 
is extremely small at A. If not for the effect of running (and in the absence of other 
contributions to i?(A)), one would find a weak scale axion, in confiict with observation. 
However, we find that a large enough B is generated in the process of running so that 
m\ = -25/i/sin2/5 ~ (400 GeV)^. 

After including the effects of renomalization group evolution and D-terms, the 
mass spectrum of the model looks as follows. The gaugino masses obey a simple 
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scaling with the the gauge couphngs, so that 



M, = M,(A) 



a,; 



(5) 



The masses of Er, II and Hd are given by 




-Ml [--si,] COS2/5 



m\^{k) + -Ml 
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where the last terms are the D-term contributions, and should be omitted for 
The A and B parameters are given by 



A~Ai(A) + M2(A) 



5~5(A)--A(A) + M2(A) 



-1.85 + 0.94 



Y 

-0.12 + 0.47-^ 
+f. 



which hold for low to intermediate values of tan (5. Here Yt = hf is the top Yukawa 
coupling squared, and Yf ~ 2.79. Au^c are obtained from At by setting Yf = 0. The 
final mass ratios for q : g : ii : W : in '■ Ni are given by 9.3 : 6.4 : 2.6 : 1.9 : 1.35 : 1.0. 
(All of our numerical results are derived using the complete one-loop RGB's.) 

In the minimal messenger model without any modification of the Higgs sector, 
one expects A and B parameters to be negligible at A, since they arise only from 
higher loops (see however Ref. 0). However, in the process of running, significant A 
and B can be induced. Furthermore, in this constrained version, the parameter tan/? 
and the sign of n are determined. These results follow from minimizing the Higgs 
potential, which leads to 



sm2/3 



fx 



-25/i 



m 



+ 2/i2 

tan2/5 - 1 



-Ml 
2 ^ 



(9) 



The last relation implies that /i is large in the model. If -B(A) = 0, we can infer its 
effective value at Mz- Tan/? is then determined through the RGE-improved tree- level 
potential of Eq. (||) to lie at tan/5 ~ 30; estimates of the full 1-loop contributions 
appear to push tan f3 up to values near tan f] ^ 40 to 50. In this most economical 
version, one also finds the sign of /i to be negative. 



5 



A further constraint arises from demanding that the mass-squared not turn 
negative (otherwise the vacuum will be superconducting). In Figure |l] we show the 
contour of = in the tan/5 — M2 plane; tan /3 ^ 50 is excluded. A more stringent 
constraint along these lines is to require, in the spirit of the phenomenology we are 
attempting to address, that > rn^^, otherwise A^^i will decay dominantly into 

T + tr, destroying the 77 signature. The contour where tr mass equals that of A^i 
is also shown in Figure |l|; above that line the tr is the lighter. Because the bound 
corresponds roughly to tan/? ~ 25, we shall focus only on the range tan/5 < 25 in 
presenting our numerical results. 

3 Production and Decay of Sparticles 

Given the physical parameters M2, tan/5 and sgn^u all cross-sections and branching 
ratios within the MMM can be calculated. It is particularly important that within 
the MMM there are three widely separated mass scales: light particles with only U{1) 
hypercharge interactions, intermediate mass particles with SU{2)l interactions, and 
heavy particles with SU{3)c interactions. The heavy states (squarks and gluinos) are 
typically around 1 TeV for light particle masses consistent with LEP bounds; there- 
fore their production rates even at the upgraded FNAL will be orders of magnitude 
below observable signal levels. Finally, the electroweak breaking constraint usually 
fixes /i such that M2 -C ^ M3, pushing sparticles whose masses come from /i 
(i.e., higgsinos and the second doublet of Higgs bosons) out of the reach of current 
accelerators. 

On the other hand, the light and intermediate mass sparticles can have non- 
negligible production rates even in Run I of FNAL (/£ ~ 120 pb"^). For this study 
of the MMM, we consider a pp collider /detector with ^/s = 1.8 TeV, as at FNAL Run 
I, with 100% efficiency for tagging all relevant final states. Because the sparticles 
are always produced in pairs, and because the gluinos and squarks are too heavy 
to be produced at reasonable x, we need only consider Drell-Yan (DY) production 
processes. Each process is calculated at the parton level for s = X1X2S, then integrated 
over the proton parton distribution functions, fq/p{xi), in this case using those of the 
CTEQ3M set 0. Previous calculations of DY sparticle production were used or 
altered for the parton level processes studied here 0. The parameter space studied 
is spanned by the three parameters (M2, tan/5, sgn/x) in the ranges 50 < M2 < 400, 
1 < tan /5 < 25 and fi taking both signs. 

There are a limited number of sparticle DY processes which can be expected to 
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be observable at the current generations of hadron colliders: 



pp 




pp 



W 




where i,j = 1,2. Not surprisingly, in order for production rates of the heavier slepton 
modes {VV, El^l) to be observable, rates for the lighter sparticles would be much 
larger, already providing an unambiguous signal at a collider. Since a large signal 
has in fact not been observed we will not consider the DY production of left-handed 
sleptons further. 

One also expects (and we have explicitly verified) that the rate for double neu- 
tralino production (NiNj) is small compared to the others, since in the limit Ni = 
and N2 = there are no ZNiNj interaction terms in the MSSM. (This fact also 
implies that the calculation of Ni mass bounds from non-observation of Z — »• NiNi at 
LEP1.5 need not coincide with the kinematic limit of approximately 65 GeV; in fact 
we find the appropriate mass bound from OPAL |10| to be 10 — 15 GeV below this 
limit.) Likewise, the NiW channel is suppressed due to the lack of a W~^W~B inter- 
action. We find the production rate of NiW to be at least two orders of magnitude 
below the rates for competing processes, making it inaccessible at current colliders. 

The dominant DY production channels are then inin, N2W, and WW. In Figure H 
we have shown contours of cross-sections for the three production channels in the 
(M2,tan/3) plane, where we have summed over slepton flavors and gaugino charges. 
The contours are labelled in femtobarns (fb), and one should recall that a = lOfb 
corresponds to 1 event in 100 pb~^ of data. One notes some general features: (i) 
There is little or no dependence on sgn fi in the most important channels, the largest 
dependence on sgn fi coming in the WNi channel, which is several orders of magnitude 
below those shown; (ii) the cross-sections are only weakly dependent on tan/?; (in) at 
large M2, the slepton production is largest, while it is the gaugino production which 
dominates at small M2; in particular, when the cross-sections are near 10 fb, all three 
cross-sections are within a factor of four of each other; (iv) for M2 ^ 160 GeV, the 
production rate for gauginos would be very large even now at the Tevatron. 

The exact signature of each channel will depend on how its constituents decay. 
Consistent with the MMM and the CDF event, we take the gravitino to be the 
lightest R-odd particle, to which the next-lightest R-odd particle (the A^i) decays via 
Ni — > 7G. Then essentially all sparticles produced will eventually decay to Ni and in 
turn to 7 +^±; thus 77 is a mandatory ingredient of the final state in all MMM 
sparticle production processes. 

Within the MMM, the £r = cr^JIr and the tr can only decay directly to the 
2-body final state INi or riVi, where i = e, fi (we treat r leptons separately for later 
convenience). Since the production cross-section for each flavor of slepton is the same. 
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Figure 2: Cross-sections at the Tevatron for inin (solid), WN2 (dashed), and WW 
(dotted) production. (Slepton flavors and gaugino charges summed over.) The con- 
tours are labelled in fb; the labels at the top of the plot label both gaugino production 
channels, while those at the bottom label the slepton production. The contours for 1 
event in 100 pb~^ are in bold. 



one observes a signal r+r^77+^_|_ 1/3 of the time, and £^£"77+-^^ 2/3 of the time. 
Such a final state describes the 6677 event at CDF. 

There has been some discussion in the literature as to whether the MMM could 
describe the CDF event through iRin production given the small mass splitting be- 
tween Er and Ni 0. The authors of Ref. were concerned because one of the 
electrons in the CDF event is highly energetic {E± ~ 59 GeV) while they assumed 
rrij^ ~ 1.1 m^^, which in the mass range of interest corresponds to — Ni mass 



splitting of only about 10 GeV. The required boost for the initial pairs of sleptons 
seemed to reduce significantly the probability of interpreting this event as coming 
from MMM iji production. This may be too pessimistic given the analysis of the 
Ir-Ni mass difference in the previous Section. From that analysis, we find for the 
appropriate mass range that mj^ — ~ 35 GeV and the required boost is small 
(7 ^ 1.3). However, as pointed out in [Q, larger mass differences in turn imply smaller 
rates for Ir^r production with respect to the WN2 and WW channels. For a 35 GeV 
mass splitting, this suppression is about a factor of two relative to a 10 GeV splitting. 
Given the current paucity of data, this is completely consistent with observation. 

Because N2 and W are nearly degenerate in mass, the preferred decay of the 
heavier of the two to the lighter via W"^ — > N2W^, cannot occur. If rriy^ — m^^ > M] 



w 



then the 2-body decay W^ — > NiW^ will be allowed. In Figure ||(a) the contour along 
which the decay becomes kinematically allowed is displayed. However, the coupling 
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is suppressed by the lack of any W W^B interaction in the MSSM, so unsuppressed 
3-body decays can compete with the 2-body mode. (Calculations of all relevant 2- 
and 3-body decay widths exist in the literature The 3-body decays proceed 

through t-channel exchange of II ot V: W iu^Ni. (In this way it might be 
possible to account for the CDF event using WW production, without having to 
invoke coincidence against a | x | suppression in the W decay branching ratios. A 
signal for this possibility is the existence of equal numbers of like and unlike lepton 
pair events, e.g., both e"'"e~77 and e~^fj,~'yj +^±.) 

At large tan/3 there is a strong enhancement in the 3-body decay into tv^-Ni. 
This arises because at large tan/3 the r- Yukawa coupling becomes 0{1) and couples 
to the Hd component of W , which can be relatively large. For example, for positive 
(negative) and tan/? greater than 5 (10), the branching ratio oiW tVt-Ni exceeds 
50%. _ 

The decay channels of N2 are similarly a mix of 2-body decays suppressed by small 
couplings ^nd 3-body decays suppressed by phase space. The three 2-body decays are 
to Nih^, N\Z and -^^-^^ final states. The first two can be either allowed or disllowed 
in alternative large regions of the parameter space; the third is always allowed. Each 
2-body process suffers from small couplings which go to zero in the limit A^2 = W^ 
and A^^i = . For small gaugino-higgsino mixing the leading 2-body coupling is 
through the r- Yukawa coupling so that even at relatively small tan/?, decays to r's 
predominate. 

When it is kinematically allowed, the decay — > has a large (~ 50%) 

branching fraction. Such a process could provide an interesting route toward pursuing 
Higgs physics at the Tevatron, with Higgs production neatly tagged by the hard 
photon pair. Because the lightest Higgs boson, receives large radiative corrections 
to its mass, there is considerable uncertainty about the precise prediction of m/^o 
for each point in our parameter space. Further, myfi is sensitive to deviations from 
minimal messenger models, for example to inclusion of a singlet Higgs field (perhaps 
useful for generating a /i-term) . Therefore in presenting our results, we have pretended 
that the Higgs decay mode is disallowed; if it is allowed at all it is appreciable. We will 
highlight where applicable the swath through parameter space where the Higgs mode 
is open in the truly-minimal MMM. When on-shell Higgs production is disallowed, the 
coupling suppression serves to make any 3-body decay through the Higgs unobservably 
small. The contour along which the Higgs decay mode becomes accessible is shown 
in Figure 0(a). 

The two 3-body decays of N2 are via t-channel £l and V: N2 i^i~Ni and 
A^2 vvNi. At large tan/3 the ttNi final state receives substantial enhancement 
from the Ho component of N2 coupling to t^t via the large r- Yukawa coupling, as 
discussed in the case of W above. This means that for large tan /3, the production 
of 1^"'"A^2 "will give a signal of r+r+r^77 nearly 100% of the time, while WW 
production will look the same with one fewer final state r-lepton. 

In Figure |^ we have folded the branching ratios into the production cross-sections 
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to show the expected production rates for various sets of "final" states in the detector. 
Here we have defined on-shell W^, Z, and to be final states; the reader can 
identify the observed final states associated with production of these particles simply 
by multiplying by their well-known branching ratios. By keeping the gauge bosons in 
the final state, one also learns about the invariant masses of their decays products. 
Specifically, Figure ^ displays the cx-BR of nine final states for M2 = 200 GeV and 
varying tan/3. The two pairs of plots show the results for both signs of /x. 

Figure ^(a) shows the five most relevant channels at small tan (3 and /i > (all 
channels that follow include 77+^^): ££ (solid), WM (dotted), Wtt (short-dashed), 
WW (long-dashed), and WZ (dot-dashed). Figure ^(c) is the same for /i < 0. 
Figure |^(b) shows the four most relevant channels at large tan P and /i > (not 
including ii): ttt (solid), rr (dotted), rii (short-dashed), and Wr (long-dashed). 
Figure ^(d) shows the same five curves for fi < 0. In all four plots, the shaded region 
is the portion of parameter space that would allow the decay N2 Nih^. 

From this analysis, it clearly emerges that if one is to interpret the CDF event 
to suggest a x BR{pp — > iniR) ~ (5 — 10) fb, then the MMM implies that there are 
several other channels which should also be accessible at the Tevatron. If tan/5 < 5 
or 10 (depending on sgn/x), then the largest available channel is Wii (for i = e,fi) 
with W decaying into 2 jets (2/3 of the time) or a single lepton (1/3 of the time): 
pp — > ^^^~jj77+-^± or pp — * £+£~£='=77-|-^_|_. Second to this rate will be WW'yj+^± 
channels (with the usual W decay branching ratios) and Wtt'J'j +^±. There is also a 
small window around 3 ^ tan (3^5 and /i > in which WZ'yy is an important 
final state. 

For large tan/3, r production of all kinds dominates. If tan/3 > 5(10), then the 
Tevatron should be producing {nT)X'~f'~f at rates comparable (or larger than) 

ee77 for n = 1, 2, 3. 

4 Discussion 

We have described a minimal implementation of a low-energy messenger model for 
supersymmetry breaking. Outstanding features of the model include: generic pro- 
duction of final states including 77 and missing energy; clean separation of scales in 
the R-parity odd spectrum, with the squarks and gluinos much heavier than the elec- 
troweak gauginos and sleptons; and very small mixing of the SU (2) x f/(l) neutralinos 
with Higgsinos and with each other. Many specific features are predicted in terms 
of very few parameters, so that the model could be rigorously tested, or ehminated, 
with quite a small amount of relevant data. If the CDF e"'"e~77 event represents a 
first indication, such data is not far beyond current reach. 

Of course in insisting upon the absolutely minimal model we may have over- 
reached. We have shown that some apparent difficulties with the minimal model are 
removed upon careful treatment of coupling renormalization and proper inclusion of 
D-terms. Nevertheless, one must remain open to such possibilities as different overall 
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Figure 3: Cross-section times branching ratio for M2 = 200 GeV with tan/? varied 
along the x-axis. Both signs of /i arc shown: (a)-(b) ^ > 0, (c)-(d) fi < 0. The nine 
most relevant channels are displayed, divided between two subfigures for each sgn jj,. 
The shading is in the region for which A^2 Nih^ is allowed in the minimal model. 
See the text for a full description. 
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mass scales for the sleptons as against the gauginos, or additional contributions to 
the parameters of the Higgs sector. 

Finally let us briefly mention another possibility for physics capable of generating 
the CDF e'''e~77 event, which if nothing else might serve as a useful foil. If there 
were heavy pion-like particles whose decay modes were dominantly 11^ —>■ W + j 
and n° — >■ ZZ or Z'j or 77, their production and decay could induce many of the 
same final states as we have analyzed above, including the CDF event. The missing 
energy in this event would be carried off by neutrinos; but in other cases one might 
have 77 without significant missing energy, so that this alternative will be readily 
distinguishable. 
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